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1 Introduction 

Computation basically is digital data processing. This requires the use of a computer 
which processes the data following certain set of instructions called a programme. Ex- 
amples are, numerical data being processed by the executable version of a FORTRAN 
or C programme, text being edited by a word processor and a visual image being ren- 
*^ dered by a graphics application. 

Qh! The input, intermediate and output data are internally expressed in terms of certain 

basic units called bits. Each bit has two possible values and 1 and a string of such 
values corresponds to the binary representation of a number just like the more familiar 
q-i decimal representation. The advantage of using the binary representation in a computer 

is that, it is relatively easy to construct devices that possess two clearly distinguishable 
states that may be used to represent the bit values. Examples are high and low voltage 
^ . states of a capacitor and the two stable states of a flip-flop circuit. 

The devices used to represent bits behave essentially as classical systems, even 
though they may be inherently dependent on quantum phenomena for their operation. 
For example, a transistor used in a flip-flop circuit works on the basis of the semi- 
conducting properties of certain materials. This stems from the quantum mechanical 
energy band structure of electrons in those materials. However, because of the large 
number of electrons involved, quantum effects due to them add up incoherently to 
produce say, a current or voltage that behaves classically. The flip-flop or the capacitor, 
consequently exists in one of the two possible stable states and not in any arbitrary 
mixture of them. Thus at any time, a processor using such devices to represent and 
store bits, can only process a particular set of data. In order to process several sets of 
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data concurrently, one has to use several such processors and run parallel data channels 
through them. This is ordinary parallel processing. 

There are on the other hand intrinsically quantum systems which have two orthog- 
onal basis states that may be used to represent the two values of a bit. Examples are 
the up and down states of a spin half object, the two orthogonal polarization states of 
a photon and two non-degenerate energy eigenstates of an atom. Let us represent the 
two states corresponding to bit values and 1 by |0) and |1) respectively. However 
these are not the only possible states for a quantum mechanical bit or qubit as they 
are called. Because of the superposition principle in quantum mechanics, an arbitrary 
linear combination a |0) + j3 1 1) with complex coefficients a and j8 is also a possible 
state. A qubit in such a state, in a sense, carries the two possible bit values simultane- 
ously. A device based on such qubits possesses the potential for massive parallelism 
that may be harnessed to construct quantum computers which are immensely more 
powerful than their classical counterparts. We expect such a computer to be particu- 
larly useful in simulating efficiently quantum systems such as an atom, which is a task 
for which classical computers are generally extremely inadequate. 

However, the catch is that a qubit, when measured at the end of a computation 
always collapses to one or the other basic state, yielding a value which is either or 
1. Thus, even though it may be possible for a quantum computer to carry out a large 
number of computations on different sets of data parallely, at the end of the day we 
obtain the result for just one of the sets. Notwithstanding this difficulty, it is possible 
with clever design of quantum algorithms and quantum devices to implement them, to 
use quantum computers to solve certain problems that are very hard to solve otherwise. 

2 Classical Gates 

As we discussed earlier, classical computation consists of processing or transformation 
of data represented by classical bits. The elementary units that process classical bits are 
called gates. Processors used in modern electronic computers use tens and hundreds of 
millions such gates. The design is modular. So we don't have to understand how these 
gates really work. It will be enough to treat them as little black boxes with specified 
inputs and the corresponding outputs. In that case we do not have to worry when the 
internal design of a gate changes, as long as the external function remains the same. 
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The classical gates are classified according to the number of inputs. 
• Single-input gates: 

(A) The NOT gate: This simply switches the value of the input bit from to 1 and 

vice-versa. 

a N0T a 

Mathematically: NOT a = 1 © a , where © indicates addition mod 2. 

(B) The FANOUT (Copy) gate: This is simply a wire carrying the input bit that 
branches out into two others carrying the same bit. 




(C) The ERASE gate: This simply erases or resets to the input bit. 



a — »-^<o 



Two-input gates: 



(A) The AND gate: 



a- 
b- 




a and b = ab 



This produces a single output bit from two input bits. The output is 1 only if both 
inputs are 1 and is otherwise. Mathematically aANDb = ab. 
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(B) The OR gate: 




aoRb 

The output in this case is only if both the inputs are and is 1 otherwise. 
(C) The XOR gate: 




a xor b 



Mathematically aXORb = a@b . The output in this case is 1 only if one of the inputs 
is 1 and the other is . Otherwise the output is . 

The two-input gates are conveniently described in terms of truth tables which dis- 
play the outputs for various possible inputs. The truth tables of AND, OR and XOR 
gates are shown together below. 
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From the table we easily verify 



aORb = {aANDb)XOR{aXORb) 
Thus the OR gate may be constructed by combining the AND and the XOR gates. 
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Next we may combine the NOT gate with either the AND gate or the OR gate to 
obtain two more gates. 

(D) The NAND gate: 



It turns out that any classical computation can be implemented by a circuit con- 
structed out of the above set of single-input and two-input gates. More surprisingly, it 
happens that we do not even need all the above gates. For example, the other gates can 
be made out of ERASE, FANOUT and NAND gates which therefore form a minimal 
universal set. 

3 Reversible Computation 

While the inputs of NOT and FANOUT gates may be reconstructed from the outputs, 
the same is not true for other classical gates. It is in this sense these other gates are not 
reversible. For each such gate the output contains one less bit than the input. Since a 
classical bit has two possible states, the phase space is reduced and there is a decrease 
in entropy which according to the Boltzmann relation equals &ln2. According to the 
second law of thermodynamics, this must be over-compensated by a corresponding 
increase in the entropy of the surrounding. This is equivalent to heat released to the 
environment at temperature T , given by 



This is exactly like the heat given off when one molecule of an ideal gas is 
isothermally compressed to half the original volume. Thus classical computation 
using irreversible gates inevitably generates heat. For present day computers this heat 
is of course negligible compared to the heat generated by other dissipative processes 
such as the flow of current through a resistance. As we will see later, one important 




(E) The NOR gate: 




AQ > kT In 2 
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way in which the corresponding quantum gates and computers based on them are 
different, is that they are always reversible. 

It is also possible to construct classical gates that are reversible. The general idea is 
to copy some of the input bits to the output so that the input bits may be reconstructed 
out of the result and the extra output bits. An elegant implementation of this idea is 

The Toffoli gate: 



► 1 


> >• — 


► 1 

— (- 


I > — 

) — 



where a = a , b' = b and c = c © ab. 

Since a = a , b = b' and c = c © ab' , the input bits may be reconstructed from the 
output bits simply by running the gate in the reverse. Thus the Toffoli gate is 
reversible. 

For b = c = 0,c'=0 and the Toffoli gate acts as a reversible ERASE gate. For 

b = 1 and c = 0, c = a and it simulates a reversible FANOUT gate. For b = 1 and 

c = \,c = \®a = NOT a and it is equivalent to a reversible NOT gate. For b = 1, 

c =c@a = a XOR c and it acts as a reversible XOR gate. Finally for c = 0, 

c = ab = a AND b and it simulates a reversible AND. Thus the Toffoli gate may be 

used to construct the reversible versions of all basic one input and two input classical 

gates. 

Most reversible gates use extra ancilla bits (i.e. auxiliary bits that are set to stan- 
dard states) at the input and produce extra garbage bits (i.e. bits other than those 
containing the results of computation) at the output. These are in fact necessary to 
ensure reversibility of those gates. However the garbage bits have to be erased or reset 
at the end of the computation (or earlier) because of the limitations on the available 
memory. One may worry that such erasure would spoil the reversibility and generate 
heat. Fortunately there is a way to erase the garbage without destroying reversibility. 

To see how it works, let us consider a somewhat simple situation. Suppose that the 
computation starts with the ancilla bits in the state and the target bit in the state x. 
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If some of the ancilla bits are required to be in the state 1 , we can always arrange that 
using NOT gates on those. Now a reversible computation produces the result r(x) and 
some garbage g(x) . So it looks like 

(x,0,0) -> (x,r(x),g(x)) 

We may remove the last garbage bit reversibly by simply reversing the computation 
(uncomputation) . But that removes the result r{x) too. So we carry an extra ancilla 
bit again in the state and reversibly copy the result of the computation to that. So it 
looks like 

(x, 0,0,0) -> (x,r(x),g{x),0) -> (x,r(x),g(x),r(x)) 

Now we run the computation proper backwards. This, of course, does not affect the 
fourth bit. So this final step is 

{x,r(x),g(x),r(x)) -> (x,0,0,r(x)) 

We have succeeded in removing the garbage reversibly without destroying the result! 

Problem 1: What reversible gate is simulated by the Toffoli gate for c = 1? 



Problem 2: Construct the truth table for the Toffoli gate i.e. a table displaying the 
values of the output bits a , b' and c for all possible values of the input bits a, b and c. 



Problem 3: The Fredkin gate is another 3-input 3-output gate just like the Toffoli 
gate. For this gate c = c. If c = then a = a and b=b i.e. nothing happens. On the 
other hand if c = 1 then a = b and b = a i.e. the two bits are swapped. Show that the 
Fredkin gate is reversible and explain how it may be used to simulate a reversible 
AND gate. 
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4 Quantum Gates 



The quantum gates transform qubits just like the way classical bits are changed by the 
classical gates. However, this involves the time evolution of a quantum system and 
according to the laws of quantum mechanics this is described by a unitary operator. 
Thus to every quantum gate corresponds a unitary operator U. So a quantum gate acts 
on an arbitrary multi-qubit state | !//;•„) as 

Win) \Wout) =U\\l/ in ) 

The input state may be reconstructed from the output by 

W in ) = U f \yf 0Ut ) 

Thus the quantum gates are always reversible. 

Because of linearity of the unitary operators, a quantum gate is described com- 
pletely by its action on a convenient basis. For a single qubit the most convenient 
choice is the one that corresponds to the possible results of measurement and 1 i.e. 
the states |0) and 1 1) . This is known as the computational basis. For multi-qubit states 
the computational basis is obtained by tensoring the single-qubit computational basis 
vectors. 

Note that, unlike the unitary evolution the measurement process leads to collapse 
of the state to one of the computational basis vectors and is therefore non-unitary. 

We describe below some important single qubit and two qubit quantum gates by 
their actions in the computational basis. 



• Single qubit gates 

(A) The Quantum NOT (X) gate: 

Just like the classical NOT gate, for the quantum NOT gate we have 

|0> — > 1 1> and 1 1 > — ^ |0> 
The corresponding unitary matrix is 




Actually there are infinitely many one qubit gates corresponding to infinitely many 
2x2 unitary matrices. We give below two more important examples. 
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(B) The Z gate: 

Its action in the computational basis is given by 

|0>->|0> and 
The corresponding unitary matrix is 




(C) The Hadamard (H) gate: 



This is described by the Hadamard transformation which in the computational basis, 
is given by 

|0>--L(|0> + |1» and |1)^-L(|0)-|1)) 
The corresponding unitary matrix is 



H = 



1 



1 1 



y/2 \ 1 -1 / 
We generally denote the Hadamard gate by the symbol 



Note that the Z and H gates do not have classical analogues. 



Problem 4: Think about a quantum y/NOT gate i.e a gate that is equivalent to a NOT 
gate when applied twice in succession on a qubit ( (VNOT) 2 = NOT ). Write down 
it's action on the computational basis states |0) and 1 1). 
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• Two qubit gates 
(A) The Controlled NOT (C-NOT) gate 

This uses a control (upper) qubit and a target (lower) qubit as inputs 



|a> + 



|b>— |a0b> 



The action in the computational basis is as shown above, where a,b G {0, 1} . The 
target is unchanged if the control is off {a = 0) and is flipped (NOTed) if the control is 
on {a = 1). The C-NOT is an important quantum gate. We illustrate below its use, by 
constructing a circuit that swaps a pair of qubits in the computational basis. 



I a > 



|b> 




|b> 



l a > 



The action of the circuit in the computational basis is 

\a,b) -+\a,a®b) -^\{a®b)®a,a®b) = \b,a@b) -> \b,b@{a@b)) = \b,a) 

Thus a and b are interchanged. Here and in the following we are representing a 
two-qubit computational basis state equivalently as 

\a,b) = \a) <S> \b) = \a) \b) 

If the target is initially off (b = 0) then the C-NOT gate copies the control qubit at 
the output, \a,0) — > \a,a). 

This copying of course works for the computational basis states. What if the control 
qubit is in a general state | \jf) = a |0) + j3 1 1) , where | a | 2 + | j3 | 2 = 1? In that case 

|y)|0) = a|0,0)+/3|l,0)->a|0,0)+/3|l,l) 



10 



On the other hand, if the control qubit is faithfully copied, the final state should be 



| y) \y) = a 2 10,0) +/3 2 11, 1) + a/3 10, 1) +/3a |1,0) 



This is the same as above only if either a = or j3 = i.e. if the input qubit is in a 
computational basis state. Would it be possible to use more complicated gates and 
circuits to copy arbitrary quantum states? The answer, following from the linearity 
and unitarity of quantum gates, is no . This result goes by the name: No Cloning 
Theorem. 

Problem 5: Design a quantum circuit that will copy the Hadamard states 
|±) = -^=(|0) ± |1)) faithfully, using the quantum C-NOT gate and four Hadamard 
gates. Show that this circuit is equivalent to a quantum C-NOT gate with the control 
and target bits interchanged. 



Problem 6: Prove the No Cloning Theorem: Any quantum copier can at best copy a 
set of mutually orthogonal states and not any arbitrary unknown quantum state, in 
the following way. Consider a unitary operator U that takes the product state | x/r) \S) 
to \ y) where |y/) is the state to be copied and \S) is some standard state. If this is 
possible for two different states \\jf) = |a) and | \\f) = |j3) then show that they must be 
orthogonal. 

Another important two qubit gate is 
(B) The Function gate (f-gate): 



The action in the computational basis is shown above. This is a generalization of the 
C-NOT gate which evaluates a Boolean function of a Boolean argument 
/ : {0, 1} — ► {0, 1}. The C-NOT corresponds to f(x) = x. For y = , the action of the 
f-gate is 



X> + 




\x,0) 



!*>/(*)) 



n 



What if the control (upper) qubit is in an arbitrary state | y) = a |0) + j3 1 1)? In that 
case the result is 



Thus the result contains the values of the function / for two possible arguments 
simultaneously. As we will see, this is the key to the quantum parallelism alluded to 
earlier. 

Problem 7: Show that the quantum f-gate is unitary, where / : {0, 1} — > {0, 1} is a 
Boolean function of a Boolean argument. 

One may be tempted to think that the quantum gates are like probabilistic classical 
gates. For example, the Hadamard gate converts the state |0) to — ^ ( 1 0) + |1)), which 
upon measurement yields the values or 1 each with probability ^. This is just like 
a classical gate that produces the result or 1 each with probability i, depending on 
say, the result of a fair coin toss. However this notion is quickly dispelled by the 
observation that a second application of the Hadamard gate will change the state back 
to |0) that yields the value with certainty. This is of course impossible with any 
classical probabilistic gate. It happens with quantum gates because of interference 
between parallel channels. 

More formally the Hadamard state -t?(|0) + |1)) is a pure state described by the 



\y) |0> = a |0,0> +P |1,0) - a|0,/(0)> +/3 11,/(1)> 




in the computational basis, whereas the output of the clas- 



sical probabilistic gate is a mixed state with the density matrix ^ 




5 Bloch Sphere Representation 



Consider a general normalized one qubit state | y) = a |0) + j8 1 1). Since 
I a I 2 + I P | 2= 1> we ma y parameterize 

a = cos-, p=e i(p sm- (0 < d < n, < <p< In) 



up to an unimportant overall phase factor. 
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The angles 6 and «p may be used as the polar angle and azimuth of a point on the unit 
sphere as shown above. This is known as the Block sphere representation 1 of one 
qubit state \ \ff} ■ The computational basis states |0) and 1 1) are then represented by the 
north and south poles respectively, of the sphere. The single qubit gates correspond to 
transformations on the Block sphere. For example, the X-gate (NOT) and the Z-gate 
correspond to rotations through % about the x and z axes respectively. 

Problem 8: Single qubit density matrix: Show that the density matrix for a single 
qubit can be expressed as 

p = l -(i+r-i?) 

where / is the 2x2 unit matrix , a, are the three Pauli matrices and ~f is an arbitrary 
(radius) vector. Prove that for pure states | ~T |= 1 and for mixed states | 7* |< 1. 
Thus the pure states are represented by points on the unit sphere whereas the mixed 
states are represented by points inside it. This is precisely the description we obtained 
above for pure states. Density matrix formalism allows it to be extended to mixed 
states. 

Problem 9: Describe the action of the Hadamard gate on the Bloch sphere. 

We conclude this section with a brief discussion about the universality of quantum 

gates. It turns out that any quantum gate (and therefore circuit) may be simulated by 

a combination of single qubit gates and just the quantum C-NOT gate. There are, of 

course, infinitely many single qubit gates, but fortunately they may all be obtained by 

combining gates that correspond to rotations through arbitrary angles about the y and 

z axes in the Bloch sphere representation. So a minimal universal set consists of gates 

implementing these rotations R y (cc), R z (ft) and the quantum C-NOT gate. 

'This is the same as the Poincare sphere representation for the polarization states of a photon, with 
|0) representing say, the left circularly polarized state and 1 1) representing the right circularly polarized 
state. 
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6 Classical Computation with Quantum Computers 



Any quantum computer with gates simulating the basic classical gates can be used 
equally well for classical computation. We of course have to consider the reversible 
classical gates, because only those have quantum counterparts. We have seen that the 
Toffoli gate is a nice classical reversible gate that can be used to simulate the basic 
classical gates. So all we need is a quantum Toffoli gate. This, in terms of its action in 
the computational basis, is just the classical Toffoli gate with the input and output bits 
replaced by the corresponding states. 



|a> 
|b> 



0- 



|b> 



|c ©ab> 



Thus a quantum computer constructed in this way would be able to do anything that a 
classical computer can do, equally efficiently. If this mimicry is all that quantum 
computers were capable of, then there would not be much point in discussing them. 
As we will see they can do much more. 



Problem 10: The quantum half adder: Using the quantum Toffoli gate and the 
quantum C-NOT gate construct a quantum circuit that uses two single-qubit 
computational basis states \x) and \y) as inputs and produces the sum state \x®y) and 
the carry state \xy) at the output. Both input and output may contain additional states. 



7 Deutsch Problem 

Consider an arbitrary Boolean function of a Boolean argument / : {0, 1} — ► {0, 1}. 
There are, of course, four such functions corresponding to two possible arguments and 
two possible values. For two of them /(0) = /(l) and these are called constant. For 
the other two /(0) / f{\) and these are called balanced. Suppose we do not know 
the function, but are given a black box or Oracle which can evaluate it and tell us the 
result. How do we decide whether the function is constant or balanced? 

To solve this problem classically, we will have to use the oracle twice to know its 
values for and 1 . David Deutsch devised a quantum algorithm and the corresponding 
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circuit to solve this problem with just one call to the oracle. This in fact was the first 
quantum algorithm ever written. The idea is to use a quantum oracle that in a sense 
evaluates /(0) and /(l) simultaneously. This naturally employs a f-gate with / being 
our function. 



|x> 

9 d°> -h»- 



In the above, the input two qubit state is 

Win) 



x)-L(|0)-|l)) = -L(|x,0)-|x,l)) 



Hence the output is \\p out ) = -^(\x,f(x)) — \x, 1 ©/(*))) 
Since f(x) = or 1, this may be written as 

Wout) = (-l) /W -^(|x,0) - \x, 1)) = (-1)/W \x) -L(|0> - |1» 

Thus the net effect is to change the state of the top qubit according to 

W-»(-i) /w \ x ) 

i.e. the value of the function gets kicked back to the phase of the state \x). 
The actual circuit used in the Deutsch's algorithm is the following. 



|o>- 
ID- 



H 



H 



H 



The two Hadamard gates before the f-gate simply converts the states of the input 
qubits according to |0) -> ^(|0) + |1) and |1) -»■ ^(|0) - |1)). Now that the lower 
qubit is in the right state for the phase-shift action of the f-gate, the upper qubit is 
transformed linearly to 

-L[(-1)/(0)| ) + (-!)/(!) |1)] 
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If the function is constant i.e. /(0) = /(l) , then this is ±^(|0) + |1)) and the final 
Hadamard gate the produces the state ± |0). On the other hand, if the function is 
balanced i.e. f(0) / /(l) , then the result after the f-gate is the state ±^(|0) - |1)) 
for the upper qubit. Now the last Hadamard gate produces the state ± 1 1 ) . Thus a 
single call to the quantum oracle followed by the measurement of the upper qubit in 
the computational basis, solves the problem. The speedup achieved over the classical 
algorithm in this case is just a factor of 2. However a similar quantum algorithm for 
solving a generalized problem that we are going to discuss next, would show the 
power of quantum computation. 

8 The Deutsch-Jozsa Algorithm 

This solves a generalization of the Deutsch problem. Let / : {0, \} m — > {0, 1} be a 
Boolean function of a n-bit integer argument and assume that we allow only those / 
that are either constant or yield for exactly half of the arguments and 1 for the rest. 
In the latter case the function is called balanced. Given an oracle that evaluates the 
function for a given argument and returns the value, the problem is to decide whether 
it is constant or balanced. 

There are of course 2" possible arguments corresponding to that many different 
n-bit integers and to solve the problem classically, we will have to get the function 
evaluated for 

-2"+l=2"- 1 + l 
2 

arguments in the worst case. This is because, with any order of evaluation, the oracle 
may return (or 1) for first half of the arguments and we would need the value of the 
function for one more argument in order to decide if it is constant or balanced. 

The computational resources required to solve the problem grows exponentially 
with the (bit) size n of the input i.e. the argument, for large n. In the standard termi- 
nology of computer science, such problems are called hard. As we will see below, the 
Deutsch-Jozsa quantum algorithm is going to make this problem very easy. 

This algorithm uses a quantum f-gate that is a generalization of the one used in the 
Deutsch algorithm. 



|y> 



|y©f(X)> 
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The action in the computational basis is identical with that for the ordinary f-gate, 
except that \X) here, is a computational basis state of a n-qubit register labelled by a 
n-bit integer X. If the bottom qubit is in the Hadamard state -1=(|0) — |1)) , then just 
as in the case of the ordinary f-gate, the state of the upper register is transformed 
according to 

|x)^(-i) /p0 \X) 

The quantum circuit used to solve the problem is the following. 



|o>. 



H ®n 













Here the upper input | O) is the computational basis state of the n-qubit register 

labelled by the n-bit zero i.e. \0) = |0) (g) |0) <g) (g) |0) {n factors) and the 

lower input is a one-qubit computational basis state |0). 

H m = H®H® <g) H {n factors) is a generalized Hadamard operator that 

applies the Hadamard transformation on each of the n factors of a n-qubit 
computational basis state of the register. 

The effect of the first generalized Hadamard gate on the input state of the register is 

nno) = ±w + ID) ^(|0) + |i» -L(|0) + |i» = ^EB 

On the other hand the Hadamard gate acting on the input state 1 1) of the lower qubit 
puts it in the state -^(|0) — |1)) which is just right for the phase-shift action of the 
f-gate. Thus the f-gate changes the state of the register to 

2"-l 



v2" x =o 



To see the effect of the final generalized Hadamard gate on this state, we need to 
know its action on a computational basis state of the register. This is given by 

H®»\X) = ^= 2 "f(-l) x - Y \Y) 



A Y=0 



17 



where X ■ Y is the bitwise scalar product defined in the following way : if 

X = x n -\ x\xq and Y = y n -\ y\ yo are the binary representations of two n-bit 

integers then 

x-Y = e'^ Xiyi 

Problem 11: Prove the above formula for H m \X). 

Thus the state of the register finally is 

4= 2 £(-i)/W//^|x) = 1 2 £ {-l)f^\Y) 
v 1 x=o z x,y=o 

Now, the amplitude of \0) in this state is ^£|"~o(-l) /(Ar) . 

If / is constant, then this is simply ±1. On the other hand if / is balanced, then one 
half of the terms in the sum precisely cancel against the other half and the result is 0. 
Hence the probability of observing O is 1 if / is constant and is if it is balanced. 

A single call to the quantum oracle followed by measurement of the register and 
checking the result for O , allows us to decide if the function is constant or balanced. 
The quantum algorithm has achieved an exponential speedup over classical computa- 
tion! 

Problem 12: Bernstein-Vazirani problem 

Given an oracle which evaluates for some n-bit integer A, the function /a(X) = A X 
of «-bit integer X, where A -X again is the bitwise scalar product, the problem is to 
determine A. If the Deutsch-Jozsa circuit is used with the function f& , then show that 
the final state of the ra-qubit register is 

yJt (-1) A ™» = |A> 

Z X,Y=0 

Thus a single use of the oracle followed by measurement of the ra-qubit register will 
yield the integer A with certainty. 
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9 Grover Search 



The problem is to search for a particular item in an unstructured or unsorted database. 
Consider, for example, the Kolkata telephone directory. It is, of course, arranged in the 
alphabetical order of names, but not in the order of telephone numbers. Thus looking 
for a particular telephone number, involves searching an unstructured database. 

It is convenient to index the database and determine the indices for matching en- 
tries (solutions). If there are N entries in the database, then they may be indexed by 
integers 0, 1,2 N — 1. Assuming that the entries occur perfectly randomly in rela- 
tion to the search field (i.e. the telephone number in our example) , the average number 
of lookups required to find a matching entry is 

11 11 N(N+\) N+l 

1.-+2.- + +N — = —.— = 

N N N N 2 2 

If for the sake of analysis, we assume that N = 2", then it scales as 2 n ~ 1 for large n. 
This grows exponentially with the (bit) size n of the database. Hence the problem is 
hard according to the standard definition. 

Grover search is a quantum algorithm that makes the search more efficient. It does 
not make it easy though. As we will see, the number of lookups required to find a 
matching entry with high probability, scales as \/N = 2 n > 2 , which is still exponential 
in n. Note that this algorithm, unlike the Deutsch and Deutsch-Jozsa algorithms, may 
not always yield the correct result. This statistical nature, in fact, is shared by many 
quantum algorithms. However, given a candidate solution, it is usually easy to check 
its correctness. For our problem, one just has to look up the entry using the solution 
index and verify that it contains the item being searched. 

We start by defining a search function / : {0, l} 0n — ► {0, 1} such that 



1 if xis a solution 
otherwise 



The search function is evaluated by a black box or oracle which need not know the 
solutions beforehand. Given an index X as the argument, it just has to look up the 
corresponding entry in the database and check if that contains the search item. 

The oracle used in the Grover algorithm is actually a quantum oracle that uses a 
generalized f-gate just as in the case of the Deutsch-Jozsa algorithm, with the function 
/ being the search function. 



19 



PO- 



|y©f(X)> 



Its action in the computational basis is shown above. The upper input labelled by the 
index X , is a computational basis state of a n-qubit index register and the lower input 
is a computational basis state of an oracle qubit. The oracle qubit is flipped if and 
only if, the index X is a solution to the search problem. 

The oracle qubit is actually initialized in the Hadamard state 
^(|0) — 1 1)). In that case the net effect is to change the state of the index register, as 
in the Deutsch-Jozsa algorithm, according to 

\X)^0\X) = (-l)fW\X) 

where O is the oracle operator. Thus the oracle marks the solution(s) by shifting the 
phase of the index state. 

The complete circuit used in the Grover algorithm is shown below. 



|o>- 



1> 



H 



®n 



H 



H)n|-|CPSH|-|® n 



The oracle qubit starts in the state 1 1), which is then changed by the Hadamard gate to 
the Hadamard state ^(|0) — |1)). This is just right for the oracle operation. Similarly 
the index register starts in the state | O) which is then transformed by the generalized 
Hadamard gate H m to the uniform state 



N-\ 



\ w)=H ®n\ 0)= y\X) 



which has the same amplitude for all index states \X). 
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Now a set of four operations is applied in the following order. 



1. The oracle operation O 

2. Generalized Hadamard transformation H m 

3. Conditional phase shift (CPS): \X) -> (-l) 1+<5 *.o \X) . This changes the phase 
of all index states but | O) by — 1 and is therefore, equivalent to the action of the 
operator 2 \0) (0\ — I, where / is the identity operator. 

4. Generalized Hadamard transformation H m 

The product of the four is called the Grover operator G. 

Note that the product of the last three is 

H m {2\0) (0\ -I)H® n = 2\y) (y\ -I 
Thus the Grover operator is G = (2\y) (y/j —1)0 

We will see that each of the two factors in G is a reflection and therefore G itself is a 
rotation. 

To show this, it is convenient to define orthonormalized states which are uniform 
superpositions of the solution states (say, M in number) and the non-solution states 
(N — Mm number) separately, 



l«) = ^L"|X)and|/3) = ^l'|X), 



where the prime and the double-prime indicate sums over solution and non-solution 
states respectively. Then under the oracle operator O 

\a)^\a) but |/5)^-|/5) 
So O is a reflection about \a) in the \a) , |j8) plane. Now can be expressed as 



n-m. , m e. , . e lo , 

M = V — jv~ |0f) + V = cos 2 + sm 2 |j3) 



where sin f = a/ ^ 
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Thus | y) is a vector in the \a) , |/$) plane and it is easy to see that 2 1 1//) (y/j — / is a 
reflection about | ty/) in that plane. 




Simple geometry shows that the effect of these two successive reflections on \ 
is to rotate it counterclockwise towards |/3) by an angle in the \a) , |j3) plane .Thus 
the Grover operator changes | iff) to 

_ i , 30 . . . 30 . 

G\y) =cos — |a)+sm— |/3) 



Problem 13: Show algebraically that the Grover operator G rotates an arbitrary state 
in the \a) , |j8) plane by an angle 6 in the counterclockwise direction. 

After k iterations of the Grover sequence the state of the index register is 

. (2jfe+l)0. . . (2Jfc+l)0 lo . 

G k | y) = cos V 2 ' \a) + sin v 2 ; |J8) 

It is clear that after sufficient number of iterations this state would be closest to the 
solution space vector |j8). Then a measurement of the index register would yield a 
solution with a high probability. 

What is the optimum number of iterations required for this? We would obviously 
need 

(2*+ 1)0 % . , % 1 

~— i.e. k~ 

2 2 20 2 

But for large databases i.e. for large values of N, 



0=2 arcsin 
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Thus k^jyjjj. This of course depends on the number of solutions M. If we already 
know that there is just one solution i.e. M = 1 (this for example, is true in our 
telephone directory example), then k ~ \ y/N and we need 0(y/N) calls to the oracle 
compared to 0(N) in the classical case. Thus the quantum algorithm achieves a 
quadratic speedup over the classical one. 



Problem 14: Consider the Grover search for one item (M = 1) in a database of four 
entries^ = 4).What is the optimum number of calls to the oracle required? What 
then is the probability of finding the correct item? 

We end this section with a few comments. 

• The Grover sequence can in principle begin with the index register in any state in 
the plane spanned by \a) and |/$). We choose the uniform state \\\r) in this plane 
because, it is easily obtained from the computational basis state | O) without any 
a priori knowledge about \a) and |j3). 

• We need to know the number of solutions M to estimate the optimum number of 
Grover iterations. This may not be known in general before the search. Fortu- 
nately, there exist quantum algorithms for determining the number of solutions 
(without actually finding them) efficiently. 

• It has been shown that the Grover search is optimal in the sense that no other 
algorithm based on a quantum oracle can do better. 

10 Phase Estimation 

Suppose U is a unitary operator that acts on k-qubit states and \u) is an eigenstate with 
eigenvalue e l § : 

U\u)=e*\u) 

The problem is to get the best n-bit estimate for the phase fraction The quantum 
circuit used to solve the problem employs two new gates which we introduce in the 
following. 
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• The Controlled-U gate 



|a>— > — t — |a> 

,x>— [in- — {iTi> 

The action in the computational basis is shown above. The target state \X) remains 
unmodified if the control (upper) qubit is off (a = 0) and a unitary operator U is applied 
to it if the control qubit is on (a = 1). 

• The Quantum Fourier Transform (QFT) gate 

Let a x be a real number labelled by a n-bit integer x . The discrete Fourier transforma- 
tion is defined by 

The quantum Fourier transformation (QFT) is the quantum analogue of the above, 
where the numbers a x are replaced by the computational basis states \X) of a n-qubit 
register labelled by n-bit integers X . 

For n = 1 this is just our old friend the Hadamard transformation. 
The inverse transformation is 

2"-l 



/On ^ 

L X=0 



Problem 15: Show explicitly that QFT is a unitary transformation. 
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The gate that implements QFT is called the QFT gate and is represented by the 
symbol 





QFT 






— *— 



The inverse transformation is implemented by running the gate backwards. 

The actual quantum circuit used to solve the phase estimation problem is shown 
below. 



10) 

10) 
10) 



H 



H- 



H 



n 

2 - 1 



|Y,> 
\\) 

|u) 



The target k-qubit state \u) is processed together with n one-qubit states (of a regis- 
ter) which act as controls, through a series of controlled-?/ v gates. The one-qubit states 
are all initialized to |0) and then transformed by Hadamard gates to the Hadamard state 
^7j(|0) + 1 1)) before being used as inputs in the controlled-?/ 2J gates. Just to see what 
happens at the output, let us follow through the action of the first controlled-?/ gate. 

-L(|0) + |1» \u) = -L(|0> \u) + |1) |«» - -L(|0> \ U ) + |1> U 2 ° \u)) 

= -L(|0> \u) +e n °* |1) \u)) = -L(|0> |1» \u) 

So the net effect is to introduce a phase difference between the two components of 
the control state. 



-L(|o> + |i»-»-^(|o>+^|i» = | ¥b > 



where |^.) = J=(|0> |1>). 
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The effects of the other controlled- U gates are similar. So the final state of the 
control register is 

® Wj)= <g> -4(io)+^ii))=^ 2 i^i y ) 



Problem 16: Prove the above relation. 

Suppose <p/2n =X/2" where X is some n-bit integer i.e. the proper fraction has 
exactly n bits. Then the output state of the n-qubit register is 

i 2"-l 



V2'- 7=0 

This is just the result of the QFT applied to the computational basis state \X). Thus 
the inverse transformation i.e the QFT gate applied backwards to the output would 
yield \X). Hence by measuring the final output in the computational basis, we would 
obtain the phase fraction ^ exactly, with certainty ! 

What if the phase fraction has more than n bits? In that case we write 

<h X 

— = h£ 

2n 2" 

where X, as before, is some n-bit integer and < £ < ^rr- Now the output state of 
the n-qubit register in the phase estimation circuit is 

1 2"-l 



E e 2m^ e 2meY \ y \ 



'2" Y=0 

Thus the result of inverse QFT applied to this state is 

2"-l 2"-l 

_ £ e 27ti ¥ e 27lie¥ <?~ 27t '^|Z) 

2" 7=0 Z=0 



i 2"— 12"— 1 
= ~n I I e 1 ™^ e 2%UY \Z) 

2" Z=0 7=0 

The amplitude of \X) in the above is 



1 21 1 (p 27tie V - 1 

2" 2" e 2me - 1 

l g i»e(2--l) sin Ore 2 B ) 



2" sin(7T£) 
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The corresponding probability is 

P(X) 



1 sin(7T£2") 
2" sin(^£) 



Since < e2 n < j, we may bound the probability using the inequality 

2x < sin Tlx < Tlx for x G [0, -] 
Thus sin(7T£2") > 2e2" , sin(7re) < ne and hence 

2 



/ 1 2e 2" \ 4 
P(X) >(±-t±±.) = — ~ .405 



So a measurement of the final output state of the control register in the computa- 
tional basis, yields the first n bits of the phase fraction with a probability better than 
40%. Actually the probability can be made at least 1 — S for any 8e (0, 1) by using 
n + log 2 (2 + ^) qubits in the control register and rounding off the result of measure- 
ment to n bits. 



Problem 17: Suppose we initialize the target register in the phase estimation circuit 
to a linear combination £cy \ uj) of the eigenstates \uj) of the unitary operator U, 
U\uj)=e^\uj) .What do you think is obtained if the final state of the control 
register is now measured in the computational basis? Assume for the sake of 
simplicity that, each phase fraction ^ contains exactly n bits. 



The phase estimation algorithm has many interesting and powerful applications. We 
describe next, two of its important uses. 

11 Order and Factorization 

We begin with an elementary result in number theory. Let N and a be positive integers 
such that a <N and a is coprime to N i.e. gcd(a,N) = 1. Then there exists a smallest 
positive integer r <N such that a r = 1 modN. 
The integer r is called the order of a modN. 

Problem 18: For a positive integer Af , show that the set of positive integers less than 
and coprime to Af form a group under multiplication modulo N. If we denote the 
order of this group by (p(N), then (p is called the Euler (p-f unction. 
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Problem 19: Show that the order r of a mod N is a factor of <p{N). Hence prove 
a <P( N ) = i modN if a < N is coprime to N. If p is a prime, then (p(p) = p — 1 and we 
get aP~ x = 1 mo<i for a < p. Show that this is also true when a > p is not a multiple 
of p. So the general result is that a p ~ l = 1 mod p if a is coprime to the prime p. This 
is known as Fermat's little theorem. Its generalization for arbitrary positive N is due 
to Euler. 



Problem 20: Find <p(28). What is the order of 5 mod 28? 

Finding the order in general, is a a hard problem in classical computation. We 
describe below a special case of the phase estimation algorithm due to Peter Shor, that 
makes it easy. 

Suppose N is a m-bit positive integer. For the given positive integer a less than and 
coprime to N, we define a unitary operator U a such that for any m-qubit computational 
basis state \X) 

U \X)-{ \ aXmodN ) if x <N 
j \X) otherwise 

Problem 21: Show that the operator U a defined above is unitary 

It is easy to see that for each k e {0, 1 , r — 1 } the state 

1 k 
\u k ) = -=Ye- l7Ci ± \a j modN) 
V r j= o 

is an eigenstate of U a with the eigenvalue e 2m r , where r <N is the order of a modN. 
Problem 22: Show that U a \u k ) = e 2ni ; \u k ) for it = 0, 1 r - 1 



We may in principle use the phase estimation circuit with a «-qubit control register to 
obtain a rc-bit estimate of the phase fraction * with a probability exceeding 40% and 
then determine the order r from that. The catch however is that, r must be known in 
order to prepare the target register in the eigenstate \u k ) . The problem is obviated by 
the observation that the uniform superposition of eigenstates 



1 



r-l 



—ft L K) = I 1 ) 

V r H) 

and this state is therefore, easily prepared. 
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If now the target register is initialized to the ra-qubit computational basis state 1 1), 
then a measurement of the output state of the ra-qubit control register would allow us 

to make a n-bit estimate of the phase fraction * for some random ke{0,\, r— 1} 

with a probability of better than 40%. 

Problem 23: Prove that L£=o g27n '~ \ u k) = modN) . Hence show that the state 
before the inverse QFT in the order finding algorithm is 

2 f\Y)uZ\\)= 2 Z\Y)\a ¥ m odN) 

Y=0 Y=0 



How do we get the order r from the ra-bit estimate of * without knowing kl The 
solution hinges on two things; the observation that * is a rational number and the 
following theorem from the theory of continued fractions. 



Definition: A convergent is a rational fraction obtained by truncating the continued 
fraction expansion of a number after a certain number of terms. 

Theorem: If * is a rational fraction and p is a positive number such that 

1 7 — p\ < jr 1 ' tnen r occurs as a convergent in the continued fraction expansion of p. 

If X is the result of measurement of the ra-qubit control register, then with better than 
40% probability, ^ is the best ra-bit estimate of In that case, it differs from * by at 
most 

1 1 1 



2«+i 2 n +^ 2" 

Thus 



k X 
r ~ 2" 



1 

< — 
- 2« 



If we choose n > 2m , then using r <N and Af < 1 m we have 

2r 2 < IN 1 < 2 2m+l < 2" 

Hence 



k X 



1 

- 27 1 



Thus the condition of the theorem is satisfied and * can be obtained in the completely 
reduced form as a convergent in the continued fraction expansion of It also turns 
out that there is a unique convergent s - such that t < r < N < 2 m and | f — | < 2 - ' 1 . 
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Once the correct convergent j is determined, then assuming that k and r are co- 
prime we have, k = s and r = t. This assumption may be incorrect. Moreover, it is 
not certain that, ^ is the best ra-bit estimate of the phase fraction. Thus the order r 
obtained in this way must always be checked for a r = 1 modN and the computation is 
repeated if necessary. 

How is the order related to integer factorization? Assume that we have found the 
order r of a modN to be even. Then 

a r = ImodN =¥■ {ah + \){ah - 1) = OmodN 

Now, ah 7^ 1 modN; otherwise the order of amodN would at most be |. Assuming fur- 
ther that ah ^ — 1 modN, we conclude that one or both of ah ± 1 share common factors 
with N which may be found by evaluating gcd(ah ± \,N) by the standard Euclidean 
algorithm. We illustrate all these with an example. 

Example: Factoring 15 

We randomly choose a number a = l, which is less than and coprime to N = 15. 
15 is a 4 bit number i.e. m = 4 and n > 2m = 8. We take n = 11 . Suppose we 
find X = 1536 to be the result of measurement of the 11-qubit control register in the 
quantum order finding circuit. Then ^ = ^| is (hopefully) the 11 -bit estimate of the 

phase fraction * for some ke {0, 1 r — 1}. We now find the continued fraction 

expansion 

1536 _ 1 
2048 "T+j 

Thus the successive convergents are 1 and |, of which | is obviously the appropriate 
convergent. Hence assuming k and r to be coprime we find r = 4. Fortunately r is even 
and we check a r = 7 4 = 1 mod 15. Also ai = 7 2 = 4 mod 15 / —1 mod 15 . Next we 
find the two factors of Af = 15 by using the Euclid's algorithm to evaluate 

gcd(ah + \,N) =gcd(50,\5) =5 and gcd(a^ - 1,N) = gcd(4S, 15) = 3 

Thus 15 = 5 x 3 ! 
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We summarize below the basic steps in the Shor factorization algorithm 

1. Given a positive integer N , we choose another positive integer a < N randomly 
and compute gcd(a,N) using the Euclidean algorithm on a classical computer. 
If gcd(a,N) / 1 , then gcd(a,N) is a factor of N and we divide ./V by it to get the 
other factor. 

2. If gcd(a,N) = 1 i.e. a is coprime to N, then we find its period r modulo N using 
the quantum order finding circuit. 

3. If r is even and ah / —ImodN, then we compute gcd{ah ± l,N) using the 
Euclidean algorithm on a classical computer to get a pair of factors of N. If not, 
then we go back and repeat the steps. 

Why is quantum factorization important? It turns out that, for large m the best clas- 

l 2 

sical algorithm for factoring a m-bit integer requires 0(e cm ' : [ ' ^ ') steps in the worst 
case, where c is a constant. This is exponentially hard. In fact the successes of many 
cryptographic algorithms, such as the famous RSA (Rivest-Shamir-Adleman) proto- 
col, depend crucially on this hardness. Shor's quantum algorithm, on the other hand 
achieves the same goal in 0(m 2+e ) steps, where £ is small, excluding the steps needed 
by the classical Euclidean algorithm for gcd and the continued fraction expansion. The 
latter need 0(m 3 ) steps and thus the problem becomes easy. 

12 Quantum Algorithms in General 

The phase estimation algorithm can be used to solve efficiently, many other problems 
such as, period finding, discrete logarithm, Abelian stabilizer etc efficiently using a 
quantum computer. Most general among them is the 

Hidden subgroup problem: Let K be a subgroup of a finitely generated group G 
and / be a function from G to a finite set X, which is constant and distinct on each 
coset of K. Given a black-box (oracle) which implements the unitary transformation 
U\g) \h) = \g) \h®f(g)) where |g)and \h) for g G G, h G X, are vectors in Hilbert 
spaces of appropriate dimensions and © is a suitable binary operation on X , find a set 
of generators for the hidden subgroup K. 

All the quantum algorithms discovered so far, except the Grover search and its general- 
izations, are in fact special cases of Kitaev's algorithm for the Abelian hidden subgroup 
problem. 
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13 Entanglement 



|x>-HI 
|y>-> — 



0- 



lfc y > 



In the above circuit, the output for x = y = is the two qubit state 
IA)o) = 7f (1 00) + 1 11))- This state can not be expressed as the product of two 
single-qubit states and the two qubits are said to be entangled. Note that, the two 
qubits are correlated in this state; they are both or 1. This correlation survives even 
if the two qubits are separated by a large distance without disturbing the state. This 
sort of non-local correlation is responsible for the famous EPR 
{Einstein-Podolski-Rosen) paradox. |/3oo) and three other similarly entangled and 
correlated two qubit states |/3 i> = ^(|01) + |10)), |/3 10 > = ^(|00) - |11)) and 
|/3n) = ^(|01) — 1 10)) appear in John Bell's analysis of the EPR paradox and are 
therefore called Bell or EPR states. 



Problem 24: Show that the four Bell states form an orthonormal basis in the Hilbert 
space of two qubits. 



Problem 25: Show, for a two qubit state which is the product of two single-qubit 
states, that the reduced density matrix for each qubit in the pair corresponds to a pure 
state. Calculate the reduced density matrix of the first qubit when two qubits are in 
the Bell state \f} xy ) . Does it represent a pure state? 



Far from being a nuisance, entanglement is actually an useful resource. We describe 
next a couple of its applications in quantum communication. 

14 Super-dense Coding 

Suppose Alice in Amsterdam would like to send two bits of information to her friend 
Bob in Boston. 2 Classically of course two separate bits have to be sent. Can it be done 
quantum mechanically by sending just a single qubit? The answer is yes, provided 

2 For alphabetical reasons the sender in this sort of scenario is always named Alice, and the receiver is 
called Bob. 
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Alice and Bob share two qubits in a Bell state, say |j3oo). If the two bits to be sent 
are 00 then Alice simply sends her qubit to Bob who now has the two qubits in the 
state |/3oo). If the two bits to be sent are 01 then Alice applies an X-gate (Quantum 
NOT) to her qubit (assumed to be the first member of the pair) and sends it to Bob who 
would now have the pair in the Bell state |j8oi). Alice similarly applies appropriate 
transformations to her qubit if the bits to be sent are the other combinations 10 or 11 
and sends it to Bob. The general result is that the two classical bits xy are coded by a 
single Bell state \p xy ) which Bob now has. This sort of coding of a number of classical 
bits by a single entangled quantum state is known as super-dense coding. 

Since the four Bell states are orthogonal, they are certainly distinguishable by an 
appropriate measurement (not necessarily in the computational basis). Bob can there- 
fore "decode" the two qubit state which he has and get two classical bits of information. 

Problem 26: What transformation(s) should Alice apply to her qubit in order to send 
the bit pairs 10 and 11? 



Super-dense coding in also useful in the detection and correction of errors in quantum 
computation. 

15 Quantum Teleportation 

Alice in Amsterdam wants to send (teleport) an arbitrary one-qubit state 
\\jf) = a\0) + /3 |1) to Bob in Boston. She can not determine the state and send the 
information to Bob for its reconstruction. That would require measurements on in- 
finitely many copies of | y) and an arbitrary unknown state of course can not be cloned 
(No Cloning Theorem ). What can she do then? 

Assume again that Alice and Bob share an entangled EPR pair of qubits in one of 
the Bell states, say |j8oo) = ^75 (1 00) + |H))- The qubit to be teleported together with 
the EPR pair starts in the three-qubit state 

Ivo) = lv) IA»> = ^No) (|oo) + |ii» +p |i> (|oo) + |ii))] 

where, by convention, the first two qubits are with Alice and the third one is with Bob. 
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Communication Channel 



Now Alice puts her two qubits through a quantum C-Not gate. This entangles the 
qubit to be teleported with Alice's part of the EPR pair. The EPR pair was of course 
entangled to begin with. So the three qubits end up in the entangled state 

| v , 1 ) = -L[a|0)(|00) + |ll))+^|l)(|10) + |01))] 

Next Alice sends the first qubit through a Hadamard gate and measures her two qubits 
in the computational basis. The state of the three qubits after the Hadamard gate is 

IVi) = ^«(|0) + |1))(|00) + |11))+/3(|0)-|1))(|10) + |01))] 

= i[|00) (« |0> +j8 |l» + 101) (« |1) +J8 |0» 
+ 110> (« |0) - J8 |1» + |11) (« |1> -/3 |0>)] 



Thus the result of Alice's measurement would be one of the pairs 00, 01, 10 or 1 1 and 
the corresponding states in which Bob's qubit would be left are 

00- |ybo> = «|0)+J8|1) 

01- |ybi> = «|1)+J8|0) 

10- |yio> = a|O>-0|l> 

11- |yii> = «|1)-J8|0) 

If Alice now communicates her result to Bob (over a classical channel such as 
telephone or e-mail or a quantum channel using super-dense coding), he would know 
how to transform the state of his qubit to | y) . If the result is 00, he does nothing 
because his qubit is already in the state 1 For other possible results, he has to apply 
an appropriate combination of X and Z gates. If the outcome of measurement is m\m2, 
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then the general result is 

Z m 'X^-\ ¥mim2 ) = \ ¥ ) 

This seems like pure quantum magic! However it has actually been achieved in 
the laboratory by teleporting a coherent photon beam, including some deliberately 
introduced noise, from one room to another. 3 

We end this section with a couple of observations. 

• Teleportation does not violate the No Cloning Theorem, as the original state | yr) 
to be teleported, is modified in the process. 

• Teleportation is consistent with the principle of special relativity, as the actual 
information is physically communicated at a speed necessarily less than that of 
light. 

Problem 27: Describe how a shared EPR pair in the Bell state |j8n) can be used to 
teleport an arbitrary single qubit state \\jf) = a |0) + j8 1 1). 

16 Measurement and Decoherence 

The measurement of a qubit involves its interaction with the measuring apparatus. This 
could, for example, result in 

|0)|>n) |0)|mo) 
1 1) \m) —* 1 1) \m\) 

where \m) is the standard state the measuring apparatus starts in and |mo) and \m\) 
respectively are its pointer states after the interaction, corresponding to the qubit being 
in the states |0) and 1 1). 4 If, however, the qubit is in an arbitrary superposition a |0) + 
j8 1 1) , then the effect of the interaction would be 

Win) = (a |0> +j3 11» \m) -► \xi/ out ) = a |0) K) +j8 11) \ mi ) 

Thus the state of the qubit gets entangled with that of the measuring apparatus. 



3 see fhttp : / / www .its.caltech. edu/~qoptics/ teleport ■ html| 
4 The transition may be due to the interaction Hamiltonian H mt = |0) (0| ®Mq + |1) (1| ®M\, where 
My and M\ are operators that act on the Hilbert space of the apparatus. 
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The qubit in this state is described by the reduced density matrix 

p = Tr M \V out ) (y 0Mf | 

= \a\ 2 |0> <0| + |/3| 2 |1) (1| + a/3* (mi | m ( >) |0) <1| + a*j3 (m | mi) |1) (0| 
/ |a| 2 aj3*(mi | mo) \ 
\ a*j3(m 1 mi) IjSI 2 J 

where 7>m indicates the (partial) trace over the states of the measuring apparatus and 
the states |mo) and |mi) are assumed to be normalized. 

The measurement actually corresponds to the observation of the apparatus in either 
of the states |mo) or |mi) and ideal discrimination would require (mi | mo) = (why?). 
In that case the reduced density matrix is 

P = ('o' 2 I/*" 2 ) 

The vanishing of the off-diagonal elements implies complete loss of coherence 5 be- 
tween the two components of the state of the qubit. In general |(mi | mo)| < 1 and 
the interaction with the measuring apparatus leads to reduction in the off-diagonal el- 
ements of the reduced density matrix, implying partial decoherence. 

A qubit is really an open system interacting with its environment. The environ- 
mental interaction decoheres the system in exactly the same way as a measuring in- 
strument. Decoherence during operation is a fundamental factor limiting the reliability 
of quantum computation. We illustrate its effect in the simple case of the Deutsch cir- 
cuit. 



Decoherence 



|0>^{H^S 



1>— KEr 



Suppose there is just a single decoherence interaction that acts on the upper (con- 
trol) qubit between the two Hadamard gates. In that case the entangled state of that 
qubit and the environment, just before the second Hadamard gate is 

^[(-i) /(0) |o)b) + (-i) /(1) |i)ki>] 

5 Note that the loss of coherence is only apparent when the qubit is considered in isolation. Coherence 
still persists in the entangled state of the qubit and the measuring apparatus. Also, decoherence does not 
explain collapse of the state vector, though the statistical property of the reduced density matrix is the 
same as that of the collapsed state. 
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where \eo) and \e\) are two normalized states of the environment. 
The second Hadamard gate transforms this state to 

Worn) = 5[(-l) /(0) (|0) + |l))k () ) + (-l/( 1 )(|0)-|l))|e 1 )] 
= 5|0)[{(-l) /(0) ko) + (-l) /(1) |e 1 )} 

+ ^|i)[{(-i) /(0) h)-(-i) /(1) ki)} 

Assuming (eo \ e\) to be real, the reduced density matrix of the qubit in this state is 6 

p = Tr env \Yout) {¥out\ 

= ^[i + (-i) /(0)+/(1) ( eo ki)]|o)(o| 

+ ^l-(-l) /(0)+/(1) ^oki)]|l)(l| 
Hence the probabilities of measuring and 1 respectively, are 

Po = ^[l + (-l) /(0)+/(1) (eoki>] 
Pi = ^[l-(-l) /(0)+/(1) ^oki)] 

If the loss of coherence is complete i.e. (eo \ e\) = , then Pq = P\ = \ (indepen- 
dent of whether / is constant or balanced) and the circuit is totally unreliable. Even if 
there is only partial decoherence, the correct result is obtained with a probability less 
than 1. Hence the computation is not reliable. 

When the environmental decoherence is a continuous process, the state of the en- 
vironment changes with time t and the overlap is typically (eo(t) \ e\ (t)) = e~ Xt . The 
time Tdecoher = j is called the decoherence time. For t » idecoher, decoherence is 
essentially complete. 

17 Devices 

Any device used for quantum computation must be able to represent the quantum infor- 
mation robustly and perform a universal set of unitary transformation corresponding to 
the basic quantum gates. Moreover, we need to prepare fiducial initial states and mea- 
sure the output in an appropriate basis. The performance of a device depends primarily 
on the ratio of the decoherence time Xdecoher to the time scale of operation of a typical 
quantum gate x op . We would of course like Xdecoher to be as large as possible. This is 
6 Tr env in the following, denotes the trace over the states of the environment. 
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achieved by reducing the coupling with the environment. However the manipulating 
devices and the measuring apparatus also couple with the system similarly. Hence re- 
ducing the coupling too much would make it more difficult to control the state of the 
device 7 and measure the result of computation. Thus a suitable compromise has to be 
worked out. We describe very briefly below, some important classes of devices that 
have been used for quantum computation with some success. 

• Optical photon devices: These use single photon sources and interferometry 
using beam splitters, phase shifters and nonlinear Kerr media for cross phase 
modulation. The qubits are represented by the spatially different states of single 
photons. 

• Cavity QED devices: These exploit the dipole coupling of single atoms to a 
few optical modes present in a high-Q cavity. The qubits are represented by two 
levels of a single atom and are manipulated using laser pulses. 

• Ion traps: These employ few ions cooled 8 and trapped using electrostatic and 
RF electromagnetic fields. The hyperfine levels of these ions and low lying 
quantized modes of vibration of the ion chain as a whole 9 are then used to rep- 
resent qubits which are manipulated by optical laser beams. 

• NMR: In this case the polarized states of nuclear spins in high magnetic fields 
are used to represent qubits which are manipulated by radio frequency pulses. 

• Quantum dots: These are microscopic boxes created inside metals, semicon- 
ductors and even small molecules that confine electrons and holes by virtue of in- 
ternal electrostatic fields. The quantized energy levels of these confined charges 
are used to store the qubits. The qubits are controlled by electrostatic gates 
(analogous to phase shifters) and single mode wave guide couplers (analogous 
to beam splitters). 



This would increase Top 

8 so as to freeze their vibrational degrees of freedom 

9 Centre of mass phonon excitations 
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